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1. Problem Description: 

The general search problem of interest in this report is one 
in which an evading target leaves datum on a random fixed course 
with constant velocity and the searcher begins searching for the 
target at a specified time late with constant velocity and at 
predetermined search points. At each search point the searcher 
stops and looks for the target. At each search point a detection 
range is used for target detection. If the target is within 
detection range of the searcher at the search point then the 
target is detected and the search is ended. If the target is not 
within detection range of the searcher at the search point then 
the target is not detected and the search continues with the 
searcher moving to the next search point. If the target is not 
detected at the last search point then the search ends. 

2. Assumptions: 

The assumptions used for computation in the Stop and Look 
Detection Algorithm are: 

a. The target leaves datum with a constant velocity at time 
= 0 . 

b. The target leaves datum on a random fixed course unifor — 
mly distributed over the range (0,2*Pi). 

c. The searcher begins searching at a specified time late 
with a constant velocity. 

d. The searcher stops and looks for the target at predeter — 
mined search points. 

e. Detection is determi ni sti c , i.e., cookie cutter. 
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3. Input Requirements: 

The coordinate system -for describing the search is a two 
dimensional (X,Y) system centered on the datum located at (0,0). 
Within this -framework the inputs required -for the Stop and Look 
Detection Algorithm computation are: 

a. The search points, 

b. The searcher velocity, 

c. The detection range, 

d. The time late, 

e. The target velocity. 

The units used for all inputs must be consistent. 

4. Modeling Procedures: 

The Stop and Look Detection Algorithm computes the 

probability of detecting the target by the searcher by evaluating 
the intersection of the target location circle and the searcher 
detection circle for each search point and totaling the 
nonoverlapping range of target course values for which a 
detection is possible. The ratio of this range to the total range 
of target course values is the probability of detection. 

When the searcher is located at search point (X,Y) at time T 
the target location circle is centered at (0,0) with radius 
T* (target speed) and the searcher detection circle is centered at 
(X,Y) with the search detection range as it's radius. 
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5. Stop and Look Detection Algorithm Description 

The Stop and Look Detection Algorithm is described as 
follows. For each search point compute: 

a. A <i > = ( Al(i> , A2(i) ). 

A<i) is the set of target course values for which detec- 
tion will occur at search point i. A<i> is determined by 
the intersection of the target location circle and the 
searcher detection circle. 

b. B(i) = UNION A < j ) for j = 1 to i . 

B(i) = C (Bli (1> ,B2i (1) > , (Bli (2) ,B2i (2) ) , (Bli (k> ,B2i (k> >3 
where 

Bli < j ) < B2i < j ) < Bli ( j + 1 ) . 

B<i) is the set of nonoverlapping target course values 
for which detection will occur at search point i. 

c. C <i ) = SUM <B2i < j ) —Bli ( j ) ) / (2*Pi ) for j = l to k. 

C<i> is the cumulative probability of detection obtained 
by the stop and look at search point i. 

The cumulative probability of detection is then computed as 
C(n). The Stop and Look Detection Algorithm computes the total 
measure of target course values for which detection will occur 
for all of the search points. Using the assumption that the 
target leaves datum on a random fixed course with constant 
velocity, the probability of detection is computed as the ratio 
of this total target course measure to 2*Pi . 
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6. BASIC Program Description: 

The IBM-PC BASIC program for computing the probability of 
detection using the Stop and Look Detection Algorithm is included 
as Appendix A. The program includes initial screen output 
information that identifies the author and provides a general 
description of the search problem. Input file and output file 
names are provided by the user and must conform to BASIC file 
name conventions. The program uses a polar coordinate system with 
the datum (0,0) as the center. Computation for the intersection 
of the target location circle and the searcher detection circle 
is made in polar coordinates. 

Inputs for the program are provided by menu controlled 
responses which will either create a new input file, use an 
existing input file or change and use an existing input file. 

□utputs from the program are provided on the screen and 
include the date, program execution time, the input values, and 
for each search point the radius of the target location circle, 
A(i), B(i) and C(i). If B(i)=B(i — 1) then B(i) is not included in 
the output. 

When there are search points for which the target location 
circle is contained entirely within the searcher detection circle 
then the program output identifies these search points by displa- 
ying A ( i ) = (0 , 2*Pi ) . In this case the probability of detection is 
1 . 0 . 

The output is also provided in an ASCII formatted output 
file. This output file may be printed on any standard printer. 
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The Stop and Look Detection Algorithm program, sample input 
and output files reside as ASCII formatted files on disk on the 
Naval Postgraduate School mainframe IBM— 3033 computer. For users 
with access to this computer these files may be transferee!. 
Information required for this transfer is: 

User Identification Number: 0105P 
Read Only Password: Andrus 

File Name File Type 



ST0PL00K BASIC 

ST0PL00K INPUT 

ST0PL00K OUTPUT 



7. Example 1: Expanding Spiral Search: 

This example computes the probability of detection for the 
searcher using an expanded spiral search procedure. The searcher 
begins the search at datum at time late with successive search 
points located around the datum in an expanded spiral. The input 
values for this example are: 



Searcher Speed: 


200 


mi les 


per 


hour 


Target Speed: 


30 


mi les 


per 


hour 


Time late: 


2 


hours 






Detection Range: 


55 


miles 






Number of search points: 
Search points: 


21 








i ( X , Y ) i 


( 


X,Y 


) 


i 



i 


< o. 


0) 


8 


( 0 , 


100) 


15 


( 0 , 


200) 


2 


( 100, 


0) 


9 


( 100, 


-100) 


16 


(-100, 


200) 


3 


( 100, 


100) 


10 


( 200, 


-100) 


17 


(-200, 


200) 


4 


< 0 , 


100) 


11 


( 200, 


0) 


18 


(-200, 


100) 


5 


(-100, 


100) 


12 


( 200, 


100) 


19 


(-200, 


0) 


6 


(-100, 


0) 


13 


( 200, 


200) 


20 


(-200, 


-100) 


7 


(-100, 


-100) 


14 


( 100, 


200) 


21 


(-200, 


-200) 
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Table 1 contains the Stop and Look Detection Algorithm 
BASIC program output -for this example. 

Figure 1.1 contains an area coverage plot o-f this example 
data including the associated search point searcher detection 
circle, the intersecting arc o-f the target location circle, and 
the gaps in coverage o-f the target course values. 

Figure 1.2 contains a plot o-f the cumulative probability of 
detection as a function of search time for this example. 

From Table 1 and Figure 1.1 it can be seen that the search 
effort at search points 1, 8, 15, 16, 18—21 is wasted and the 
effort at search points 10, 11, 13, 17 is duplicated. 

B. Example 2: Expanding Spiral Search: 

This example uses the same data as Example 1 with the 
searcher speed increased to 290 miles per hour. 

Table 2 contains the Stop and Look Detection Algorithm 
BASIC program output for this example. 

Figure 2.1 contains an area coverage plot of this example 
data including the associated search point searcher detection 
circle, the intersecting arc of the target location circle, and 
the gaps in coverage of the target course values. 

Figure 2.2 contains a plot of the cumulative probability of 
detection as a function of search time for this example. 

From Table 2 and Figure 2.1 it can be seen that the search 
effort at search points 1, 10, 13 is wasted and the effort at 

search points 11,12, 15-21 is duplicated. 



6 



9. Example 3: Grid Search: 

This example computes the probability o-f detection for the 
searcher using a grid search procedure. The input values for this 
example are: 



Searcher Speed: 


300 


miles 


per 


hour 


Target Speed: 


30 


mi les 


per 


hour 


Time late: 


2 


hours 






Detection Range: 


55 


mi les 






Number of search points: 
Search points: 


20 









i 


( X , Y > 


i 


( X, 


Y ) 


i 


( X, 


Y ) 


i 


(-100, 200) 


8 


( 0, 


0) 


15 


( 100, 


-200) 


2 


(-100, 100) 


9 


< o. 


100) 


16 


( 200, 


-200) 


3 


(-100, 0> 


10 


( 0, 


200) 


17 


( 200, 


-100) 


4 


(-100,-100) 


11 


( 100, 


200) 


18 


( 200, 


0) 


5 


(-100,-200) 


12 


( 100, 


100) 


19 


( 200, 


100) 


6 


( 0,-200) 


13 


( 100, 


0) 


20 


( 200, 


200) 


7 


( 0,-100) 


14 


( 100, 


-100) 








Table 


3 contains 


the Stop and 


Look 


Detection Algor i 



BASIC program output for this example. 

Figure 3.1 contains an area coverage plot of this example 
data including the associated search point searcher detection 
circle, the intersecting arc of the target location circle, and 
the gaps in coverage of the target course values. 

Figure 3.2 contains a plot of the cumulative probability of 
detection as a function of search time for this example. 

From Table 3 and Figure 3.1 it can be seen that the search 
effort at search points 1, 2, 5, 6, 8, 11, 13, 16 is wasted and 
the effort at search point 20 is duplicated. 
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Andrus Stop-Look Detection Algorithm 



05-24-1985 



NUMBER OF (X,Y) = 


21 


I 


TIME LATE 


2.00 


X i Y 


SEARCH VELOCITY = 


200.00 


TIME 


TARGET VELOCITY = 


30.00 


TLCR 


DETECTION RANGE = 


55.00 


A1 ,A2 



C 



SEARCH POINT NUMBER 
SEARCH POINT COORDINATES 
TIME SEARCHER AT (X,Y) 

TARGET LOCATION CIRCLE RADIUS 
TARGET/SEARCH CIRCLE INTERSECTION ANGLES 
CUMULATIVE PROBABILITY OF DETECTION 



I 




X 


Y 


TIME TLCR A 1 


A2 


C 


1 




0.00 


0.00 


2.00 60.00 0.00 


0.00 


0.00 


2 




100.00 


0.00 


2.50 75.00 5.71 


0.57 


0. 18 


3 




100.00 


100.00 


3.00 90.00 0.61 


0.96 


0.24 


4 




0.00 


100.00 


3.50 105.00 1.03 


2.11 


0. 41 


5 




100.00 


100.00 


4.00 120.00 1.96 


2.75 


0.51 


6 




100.00 


0.00 


4.50 135.00 2.77 


3.51 


0 . 63 


7 




•100.00 


-100.00 


5.00 150.00 3.55 


4.30 


0.75 


8 




0.00 


-100.00 


5.50 165.00 0.00 


0.00 


0.75 


9 




100.00 


-100.00 


6.00 180.00 5.25 


5.74 


0.82 


10 




200.00 


-100.00 


6.50 195.00 5.59 


6. 04 


0.82 


11 




200.00 


0.00 


7.00 210.00 6.02 


0.26 


0.82 


12 




200.00 


100.00 


7.50 225.00 0.22 


0.71 


0.83 


13 




200.00 


200.00 


8.00 240.00 0.65 


0.92 


0.83 


14 




100.00 


200.00 


8.50 255.00 0.92 


1.30 


0.84 


15 




0.00 


200.00 


9.00 270.00 0.00 


0.00 


0.84 


16 




•100.00 


200.00 


9.50 285.00 0.00 


0.00 


0.84 


17 




•200.00 


200.00 


10.00 300.00 2.18 


2.54 


0.84 


18 




■200.00 


100.00 


10.50 315.00 0.00 


0.00 


0.84 


19 




•200.00 


0.00 


11.00 330.00 0.00 


0.00 


0.84 


20 




•200.00 


-100.00 


11.50 345.00 0.00 


0.00 


0.84 


21 




•200.00 


-200.00 


12.00 360.00 0.00 


0.00 


0.84 


B ( 


I 


) = CUMULATIVE ANGLE DETECTION COVERAGE AT 


X (I) ,Y(I) . 




B ( 


2 


)= (0. 


00,0.57) 


(5.71,6.28) 






B ( 


3 


)= (0. 


00,0.57) 


(0.61,0.96) (5.71,6.28) 






B ( 


4 


)= (0. 


00,0.57) 


(0.61,0.96) (1.03,2.11) 


(5.71,6.28) 




B( 


5 


) = (0. 


00,0.57) 


(0.61,0.96) (1.03,2.75) 


(5.71,6.28) 




B ( 


6 


)= (0. 


00,0.57) 


(0.61,0.96) (1.03,2.75) 


(2.77,3.51) 


(5.71,6.28) 


B ( 


7 


) = (0. 


00,0.57) 


(0.61,0.96) (1.03,2.75) 


(2.77,3.51) 


(3.55,4.30) 






(5. 


71,6.28) 








B ( 


9 


)= (0. 


00,0.57) 


(0.61,0.96) (1.03,2.75) 


(2.77,3.51) 


(3.55,4.30) 






(5. 


25,6.28) 








B ( 


12 )= (0. 


00,0.96) 


(1.03,2.75) (2.77,3.51) 


(3.55,4.30) 


(5.25,6.28) 


B ( 


14 )= (0. 


00,2.75) 


(2.77,3.51) (3.55,4.30) 


(5.25,6.28) 





Execution Time = 0.46 Minutes 



Example 1: BASIC Program Output - Expanding Spiral Search 



Table 1 
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Example 1: Expanding Spiral Search Area Coverage 

Fi gure 1 . 1 



9 



CUNULAT IUE PROBABILITY OF DETECTION 



1.0 



.3 
.8 
. 7 
. 6 



.2 

. 1 
0.0 J — F 



6 7 8 3 

TINE 



10 11 12 



Example 1: Cumulative Probability o-f Detection vs Time 

Figure 1.2 
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Andrus Stop-Look Detection Algorithm 



05-24-1985 



NUMBER OF (X,Y) 


= 


21 


I 


SEARCH POINT NUMBER 


TIME LATE 


= 


2.00 


>- 

X 


SEARCH POINT COORDINATES 


SEARCH VELOCITY 


= 


290.00 


TIME 


TIME SEARCHER AT (X,Y) 


TARGET VELOCITY 


= 


30.00 


TLCR 


TARGET LOCATION CIRCLE RADIUS 


DETECTION RANGE 


= 


55.00 


A 1 , A2 
C 


TARGET/SEARCH CIRCLE INTERSECTION ANGLES 
CUMULATIVE PROBABILITY OF DETECTION 



I 




X 


Y 


TIME TLCR A1 


A2 


C 


1 




0.00 


0.00 


2.00 60.00 0.00 


0.00 


0.00 


2 




100.00 


0.00 


2.34 70.34 5.72 


0.56 


0.18 


3 




100.00 


100.00 


2.69 80.69 0.00 


0.00 


0. 18 


4 




0.00 


100.00 


3.03 91.03 0.99 


2. 15 


0.36 


5 




100.00 


100.00 


3.38 101.38 2.04 


2.67 


0.45 


6 




100.00 


0.00 


3.72 111.72 2.63 


3.66 


0.60 


7 




100.00 


-100.00 


4.07 122.07 3.53 


4.32 


0.71 


8 




0.00 


-100.00 


4.41 132.41 4.32 


5. 10 


0.83 


9 




100.00 


-100.00 


4.76 142.76 5.11 


5.89 


0.93 


10 




200.00 


-100.00 


5.10 153.10 0.00 


0.00 


0.93 


11 




200.00 


0.00 


5.45 163.45 6.06 


0.23 


0.93 


12 




200.00 


100.00 


5.79 173.79 0.35 


0.58 


0. 93 


13 




200.00 


200.00 


6.14 184.14 0.00 


0.00 


0.93 


14 




100.00 


200.00 


6.48 194.48 0.88 


1.33 


0.95 


15 




0.00 


200.00 


6.83 204.83 1.30 


1.84 


0.95 


1 b 




•100.00 


200.00 


7.17 215.17 1.79 


2.28 


0.95 


17 




•200.00 


200.00 


7.52 225.52 ' 0.00 


0.00 


0.95 


18 




•200.00 


100.00 


7.86 235.86 2.44 


2.91 


0.95 


19 




•200.00 


0.00 


8.21 246.21 3.01 


3.28 


0.95 


20 




•200.00 


-100.00 


8.55 256.55 3.42 


3.79 


0.95 


21 




•200.00 


-200.00 


8.90 266.90 3.74 


4. 12 


0.95 


B ( 


I 


)= CUMULATIVE l 


ANGLE DETECTION COVERAGE AT 


X ( I ) , Y ( I ) . 




B ( 


2 


)= (0. 


00,0.56) 


(5.72,6.28) 






B< 


4 


)= (0. 


00,0.56) 


(0.99,2.15) (5.72,6.28) 






B< 


5 


)= (0. 


00,0.56) 


(0.99,2.67) (5.72,6.28) 






B( 


b 


)= (0. 


00,0.56) 


(0.99,3.66) (5.72,6.28) 






B ( 


7 


)= (0. 


00,0.56) 


(0.99,4.32) (5.72,6.28) 






B ( 


8 


)= (0. 


00,0.56) 


(0.99,4.32) (4.32,5.10) 


(5.72,6.28) 




B ( 


9 


)= (0. 


00,0.56) 


(0.99,4.32) (4.32,5.10) 


(5.11,6.28) 




B ( 


12 )= (0. 


00,0.58) 


(0.99,4.32) (4.32,5.10) 


(5.11,6.28) 




B ( 


14 )= (0. 


00,0.58) 


(0.88,4.32) (4.32,5.10) 


(5.11,6.28) 





Execution Time = 0.47 Minutes 



Example 2: BASIC Program Output - Expanding Spiral Search 

Table 2 
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Example 2: Expanding Spiral Search Area Coverage 

Figure 2. 1 
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CUMULAT IUE PROBABILITY OF DETECTION 



1.0 



. 3 
.8 
. 7 

. 6 ' 

.5 

. 4 
.3 
.2 
. 1 

0.0 



_ y 



2 3 4 5 



6 7 8 

TlflE 



3 10 11 12 



Example 2: Cumulative Probability o-f Detection vs Time 

Figure 2.2 
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NUMBER 


OF ( X , Y ) 


= 20 


I 




SEARCH 


POINT 


NUMBER 




TIME LATE 


= 2.00 


X, 


Y 


SEARCH 


POINT 


COORDINATES 


SEARCH 


VELOCITY 


= 300.00 


TIME 


TIME SEARCHER 


AT ( X , Y ) 




TARGET 


VELOCITY 


= 30.00 


TLCR 


TARGET 


LOCATION CIRCLE 


RADIUS 


DETECTION RANGE 


= 55.00 


A 1 


,A2 


TARGET/SEARCH 


CIRCLE INTERSECTION ANGLES 








C 




CUMULATIVE PROBABILITY 


OF DETECTION 


I 


X 


Y 


TIME 


TLCR 


A1 


A2 


C 


1 - 


100.00 


200.00 


2.00 


60.00 


0.00 


0.00 


0.00 


2 - 


100.00 


100.00 


2.33 


70.00 


0.00 


0.00 


0.00 


3 - 


100.00 


0.00 


2.67 


80.00 


2.56 


3.72 


0.18 


4 - 


100.00 - 


100.00 


3.00 


90.00 


3.75 


4. 10 


0.24 


5 - 


100.00 - 


200.00 


3.33 


100.00 


0.00 


0.00 


0.24 


6 


0.00 - 


200.00 


3.67 


110.00 


0.00 


0.00 


0.24 


7 


0.00 - 


100.00 


4.00 


120.00 


4.24 


5.18 


0.39 


8 


0.00 


0.00 


4.33 


130.00 


0.00 


0.00 


0.39 


9 


0.00 


100.00 


4.67 


140.00 


1.25 


1.89 


0.49 


10 


0.00 


200.00 


5.00 


150.00 


1.44 


1.70 


0.49 


11 


100.00 


200.00 


5.33 


160.00 


0.00 


0.00 


0.49 
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100.00 


100.00 


5.67 
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0.48 


1.09 


0.59 
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0.00 
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0.00 


0.00 


0.59 
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100.00 
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0.73 
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200.00 


0.00 


7.67 


230.00 


6.07 


0.22 


0.80 


19 


200.00 


100.00 


8.00 


240.00 


0.24 


0.69 


0.84 


20 


200.00 


200.00 


8.33 


250.00 


0. 62 


0.95 


0.84 
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Execution Time = 0.49 Minutes 



Example 3: BASIC 



COVERAGE AT X ( I ) ,Y(I) . 



(4.24,5. 18) 
(3.75,4. 10) 
(2.56,3.72) 


(4.24,5. 18) 
(3.75,4. 10) 


(4.24,5. 18) 


(2.56,3.72) 


(3.75,4. 10) 


(4.24,5. 18) 


(2.56,3.72) 


(3.75,4. 10) 


(4.24,5.66) 


(2.56,3.72) 


(3.75,4. 10) 


(4.24,6.07) 


(1.25,1.89) 


(2.56,3.72) 


(3.75,4. 10) 


(1.25,1.89) 


(2.56,3.72) 


(3.75,4.10) 


ogram Output 


- Grid Search 





Table 3 
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Example 3: Grid Search Area Coverage 

Figure 3. 1 
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ClinULAT [ UE PROBABILITY OF DETECTION 




Example 3: Cumulative Probability of Detection vs Time 

Figure 3.2 
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CUMULATIUE PROBABILITY OF DETECTION 



1.0 



. 3 
.8 
.7 
. 6 
.5 
.4 
.3 
.2 
.1 
0.0 

2 3 4 5 6 7 8 3 10 11 12 

TIME 




Ex 1 

Ex 2 

Ex 3 



Cumulative Probability o-f Detection vs Time -for Examples 1,2,3 

Figure 4 



10. Assumption Modifications: 

The assumptions for this version of the Stop and Look algo- 
rithm can be changed with minor modifications to the algorithm 
procedures and the program contained in Appendix A. Examples of 
these assumption changes are: 

a. Searcher velocity: Searcher velocity does not have to be 
constant. What is required is that the time at each search point 
be known. 

b. Target course: Target course does not have to range over 
(0,2*Pi) or be uniformly distributed. If target course is distri- 
buted uniformly then C(i) is computed using the assigned target 
range. If the target course values are not uniformly distributed 
then C(i) would be computed by weighting the components of B(i) 
according to the target course distribution. 

c. Target velocity: Target velocity does not have to be 
constant. What is required is that target velocity as a function 
of time is known. 

Other assumption changes can be suggested. For example it is 
possible to assume distributions for target and searcher veloci- 
ties and compute C(i) by weighting the B(i) components according 
to the velocity distributions. However, as the assumption changes 
become more complex the algorithm modifications increase in 
complexity. 
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Appendix A: BASIC Program Listing: 



10 'ANDRUS STOP AND LOOK ALGORITHM: VERSION 1.0: 26 OCTOBER 1984 
20 DEFDBL A-H.P-Z: DEFINT 1-0: KEY OFF 

30 CLS: DEF FNACOS (A) =1 . 570796-ATN (A/SQR ( 1-A*A) ) : PI =3 . 141592653589793# 

40 DIM S(51) f ST (51 > , TR (51 ) ,TI (51) ,A1 (102) , A2 ( 102) ,X (51) , Y ( 5 1 ) ,B1 (102) , 

B2 ( 102 ) , C (51 ) 

50 Al*="Stop and Look Detection A1 gor i thm" : A2*="Wr i tten by:" 

60 A3*="Prof essor Alvin Andrus, Code 55As" : A4*="Naval Postgraduate School" 



70 A 5 *="Monterey, Ca. 93940" : A6*="To Continue use ENTER key" 

80 LOCATE 10 , 40-LEN ( A 1 $ ) /2 : PRINT All: LOCATE 12, 40-LEN (A2*) /2: PRINT A2* 

90 LOCATE 14 , 40-LEN (A3*) /2: PRINT A3*: LOCATE 16, 40-LEN (A4*) /2: PRINT A4* 

100 LOCATE 1 8 , 40-LEN ( A5*) /2: PRINT A5 *: LOCATE 24, 40-LEN (A6*) /2: PRINT A 6 *; 

110 INPUT "",A*: CLS: LOCATE 3,1 

This program computes the probability of detection for a" 
searcher looking discretely for an evading target. ".‘PRINT 
Assumptions are:": PRINT 

1) The datum is located at (0,0)." 

2) Detection is deterministic, i.e., cookie cutter." 

3) The searcher begins searching at time late with constant velocity 

4) The target leaves datum on a random course with constant velocity 

5) The searcher stops and looks for the target at specified points." 



120 PRINT" 
130 PRINT" 
140 PRINT" 
150 PRINT" 
160 PRINT" 
170 PRINT" 
180 PRINT" 
190 PRINT" 
200 PRINT 
210 PRINT" 
220 PRINT" 
230 PRINT" 
240 PRINT" 
250 PRINT" 
260 PRINT" 
270 PRINT" 



Input Requirements are:":PRINT 

1 ) The time late. " 

2) The target speed . * 

3) The searcher speed." 

4) The detection range." 

5) The stop and look Points: X ( I ) , Y ( I ) . I<=50.": PRINT 
The units used for speed, time and range must be consistent." 

280 LOCATE 24 , 40-LEN (A6* ) /2 : PRINT A6*;: INPUT "",A*: CLS: LOCATE 5,1 
290 PRINT" To Create a New Input File: Type N or n." 

300 PRINT" To Use an Existing Input File: Type E or e." 

310 PRINT" To Change and Use an Existing Input File: Type C or c." 

320 A*=INKEY* 

330 IF A*=“C" OR A*="c" OR A*="E" OR A*="e" OR A*="N" OR A*="n" THEN 350 
340 GOTO 320 
350 PRINT: INPUT " 

360 PRINT: INPUT " 



Input File Name 
Output File Name 



= "j INFILE* 

= ": 0UT1*: 0UT2*=0UT1*+" . 1 " : PRINT 



370 IF A*="E" OR A*="e" OR A*="C" OR A*="c" THEN 680 
380 PRINT" If value is correct use ENTER key." 

390 PRINT" If value is incorrect enter new val ue. ": PRINT 

400 PRINT" Searcher speed = "; S V ; : INPUT SV*:IF SV*<>"" THEN SV=VAL(SV*) 

410 PRINT" Target speed =" ; TV; : INPUT TV*: IF TV*<>"" THEN TV=VAL(TV* 

420 PRINT" Time late ="; TL; : INPUT TL*: IF TL*<>"" THEN TL=VAL (TL* 

430 PRINT" Detection range = " ; R ; : INPUT Rl: IF R*<>"" THEN R=VAL (R*) 

440 PRINT" Number of stop and look points (<=50) = " ; N ; 

450 INPUT N*: IF N*<>"" THEN N=VAL (N* ) 

460 CLS: IF N>=10 THEN 530 

470 FOR 1=1 TO N 

480 LOCATE 1+2, 1: PRINT " X " ; RI GHT* ( STR* ( I ) , 1 ) ; " = " ; X ( I ) j : I NPUT X* 

490 IF X*<>"" THEN X(I)=VAL(X*) 

500 LOCATE 1 + 2,40: PRINT "Y" 5 RIGHT* (STR* ( I ) , 1 ) ; " = " 5 Y ( I ) ; : INPUT Y* 

510 IF Y * < > " " THEN Y ( I ) =VAL ( Y* ) 
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520 NEXT I: GOTO 650 



FOR 1=1 TO 9 

LOCATE 1+2, 1: 
IF XSO"" THEN 
LOCATE 1+2,40: 
Y $ < > " " THEN 



530 
540 
550 
560 

570 IF 

580 NEXT I 
590 FOR 1=10 TO N 
600 LOCATE 1+2, 1: 

IF XSO’"' THEN 
LOCATE 1+2,40: 
YS< > " " THEN 



610 
620 
630 IF 
640 NEXT I 



1 ) 

1 ) 



640 NEXT I 

650 OPEN INFILES FOR OUTPUT 
660 PRINT#1 , N; TL; SV; TV; R: F( 



PRINT " X " ; RIGHTS ( STRS ( I ) , 

X ( I ) =VAL ( XS) 

PRINT " Y" ; RIGHTS (STRS C I ) , 

Y ( I) =VAL (YS) 

PRINT "X" ;RIGHTS(STRS(I) ,2); 
X ( I ) =VAL (XS) 

PRINT " Y" ; R I GHTS ( STRS d ) ,2) ; 
Y ( I ) =VAL C YS) 



I " = '' ; X ( I ) ; : INPUT XS 
! " = ";Y(I) ; : INPUT YS 

"=" ; X ( I ) ; : INPUT XS 
" = ''; Y(I) ;: INPUT YS 



670 CLOSE# 1: GOSUB 1500 

680 OPEN INFILES FOR INPUT AS 1 

690 INPUT # 1 ,N,TL,SV,TV,R: FOR 1 = 1 TO N: I NPUT# 1 , X ( I ) , Y ( I ) : NEXT I 
700 CLOSE# 1: IF AS="C“ OR AS="c” THEN 380 ELSE GOSUB 1500 
710 BT=TIMER: S ( 1 ) =SQR ( X ( 1 ) A 2+Y ( 1 ) A 2) : ST ( 1 ) =0: TI ( 1 ) =TL: TRd)=TL*TV 
720 FOR 1=2 TO N 

730 S(I)= SQR ( X ( I ) A 2+Y d ) A 2 ) 

740 ST ( I ) =ST ( I - 1 ) +SQR ( (X (I)-X(I-l) ) A 2+ ( Y ( I ) -Y ( 1-1 ) ) A 2) 

750 TI ( I ) = ST ( I ) /SV+TL: TRd)=TId)*TV 

760 NEXT I 

770 FOR 1=1 TO N 

780 IF R< =ABS ( TR ( I ) -S ( I ) ) THEN A1(I)=0: A2(I)=0: GOTO 910 

790 IF R<S ( I ) +TR ( I ) THEN 810 

800 A1 ( I ) =0: A2 ( I ) =2*P I : GOTO 910 

810 IF X ( I ) =0 AND Y ( I ) >0 THEN A = P I / 2 : GOTO 860 

820 IF X ( I ) =0 AND Y(IKO THEN A=3*PI/2: GOTO 860 

830 A=ATN ( Y ( I ) / X ( I ) ) 

840 IF X(I)<0 THEN A=A+PI 

850 IF X ( I ) >0 AND Y(I)<0 THEN A=A+2*PI 

860 B=(R A 2-Sd) A 2-TRd) A 2)/(-2*Sd)*TRd)) : A 1 (I ) =FNACOS ( B ) 

870 A2 (I ) =2*PI-A1 (I ) +A: A1(I)=A1(I)+A 

880 IF Aid) >2#PI THEN A1 (I ) =A1 (I ) -2*PI : GOTO 880 

890 IF A2 (I ) >2#P I THEN A2 (I ) =A2 (I ) -2*PI : GOTO 890 

900 IF A< A 1 ( I ) OR A>A2 ( I ) THEN SWAP A1(I),A2(I) 

910 NEXT I 
920 K=0 

930 FOR 1=1 TO N 

940 IF Aid) =A2 ( I ) THEN 1470 

950 IF K <> 0 THEN 980 

960 IF A 1 (I X A2 (I ) THEN B 1 ( 1 ) =A1 (I ) : B2 ( 1 ) =A2 (I ) : C (I ) =A2 (I ) -A 1 (I ) : K= 1 : 

GOTO 1470 



970 B1 (1 ) =0: B2 ( 1 ) =A2 (I ) : B1 (2) =A1 (I ) : B2 (2) =2*PI : C (I ) =A2 (I ) -A1 (I ) +2*PI : K 
GOTO 1470 

980 IF A1 (I X A2 (I ) THEN A=A Id): B=A2 (I ) : FLAG= 1 : GOTO 1000 
990 A=0: B=A2 (I ) : C=A1(I): D=2*PI: FLAG=2: GOTO 1000 

1000 FOR J=1 TO K 

1010 IF A<=B1 (J) THEN IX = 1: I X 1 = J : GOTO 1050 

1020 IF A<=B2 ( J) THEN IX=2: I X 1 = J : GOTO 1050 

1030 NEXT J 

1040 B1(K+1)=A: B2(K+1)=B: C(I)=C(I)+B-A: K=K+ 1 : GOTO 1460 
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1050 FOR J = 1 TO K 

1060 IF B<= B1(J) THEN I Y= 1 : IY1=J: GOTO 1100 

1070 IF B< = B2 ( J ) THEN IY=2: IY1=J: GOTO 1100 

1080 NEXT J 
1090 IY=3 

1100 IF IX = 1 AND I Y=1 AND IX1=IY1 THEN 1200 

1110 IF IX = 1 AND I Y= 1 AND IXIOIYI THEN 1230 

1120 IF I X= 1 AND I Y=2 AND IX 1=IY1 THEN 1260 
1130 IF IX = 1 AND I Y=2 AND IXIOIYI THEN 1280 

1140 IF IX=2 AND I Y= 1 AND IX1=IY1 THEN 1320 

1150 IF IX=2 AND IY=1 AND IXIOIYI THEN 1330 

1160 IF I X = 2 AND IY-2 AND I X 1 = 1 Y1 THEN 1460 

1170 IF I X =2 AND IY=2 AND IXIOIYI THEN 1370 

1180 IF I X = 1 AND I Y = 3 THEN 1420 

1190 IF I X =2 AND IY=3 THEN 1440 

1200 C ( I ) =C ( I ) +B-A 

1210 FOR L=K TO 1X1 STEP -1: B1 (L+l ) =B1 (L) : B2 (L+l ) =B2 (L) : NEXT L 
1220 B1 ( IX 1 ) =A: B2 ( I X 1 ) =B : K=K+1: GOTO 1460 

1230 C ( I ) =C ( I ) +B-A: FOR L=IX1 TO IYl-ls C ( I ) =C < I ) +B1 (L) -B2 (L) s NEXT L 
1240 FOR L= 1 TO K-IY1: B1 ( IX1+L) =B1 ( I Yl + L) : B2 ( I X 1+L) =B2 ( I Yl+L) : NEXT L 

1250 B 1 ( I X 1 ) = A : B2 ( I X 1 ) =B: K=K- IY1 + I X 1 : GOTO 1460 

1260 C(I)=C(I)+B1 (IXl)-A 

1270 B1 ( I X 1 ) =A: GOTO 1460 

1280 C ( I ) =C ( I ) +B1 (IYl)-A: FOR L=IX1 TO IYl-ls C ( I ) =C ( I ) +B1 (L) -B2 (L) : NEXT L 
1290 B1(IX1)-A: B2 ( IX 1 ) =B2 ( IY 1 ) 

1300 FOR L = 1 TO K=IY1: B 1 ( IX 1 + L) =B1 < I Yl+L) : B2 < IX 1+L) =B2 ( I Yl+L) : NEXT L 
1310 K=K-I Y 1 + 1 X 1 : GOTO 1460 

1320 CLS: PRINT "ERROR IN ALGORITHM LOGIC: PROGRAM EXIT": END 
1330 C ( I ) =C < I ) +B-B2 < IX 1 ) 

1340 FOR L=IX 1+1 TO IYl-ls C ( I ) =C ( I ) +B 1 (L) -B2 (L) : NEXT L 

1350 FOR L=1 TO K+1-IY1: B1 (IX1+L)=B1 (IY1+L-1) : B2 ( IX 1+L) =B2 < I Yl+L-1 ) : NEXT L 
1360 B2 ( I X 1 ) =B: K-K-IY1+IX 1 + 1 : GOTO 1460 
1370 C ( I ) =C ( I ) +B1 (IY1)-B2(IX1) 

1380 FOR L=IX 1 + 1 TO IY1-1: C ( I ) =C ( I ) +B1 (L ) -B2 <L) : NEXT L 
1390 B2 ( I X 1 ) =B2 ( I Y 1 ) 

1400 FOR L = 1 TO K-IY1: B1 ( IX1+L) =B1 ( IY1+L) : B2 ( I X 1+L) = B2 ( IY1+L) : NEXT L 
1410 K=K-I Y1+IX1 : GOTO 1460 

1420 C (I) =C ( I) +B-A: FOR L= IX 1 TO K: C ( I ) =C ( I ) +B1 (L) -B2 (L) : NEXT L 
1430 B1 ( IX 1 ) =A: B2(IX1)=B: K=IX1: GOTO 1460 

1440 C ( I ) =C ( I ) +B-B2 ( I X 1 ) : FOR L=I X 1 + 1 TO K: C ( I ) =C ( I ) +B1 (L) -B2 (L) : NEXT L 
1450 B2 ( IX 1 ) =B: K= I X 1 : GOTO 1460 

1460 IF FLAG=1 THEN 1470 ELSE A=C: B=D: FLAG=1 : GOTO 1000 
1470 IF 1 = 1 THEN C ( 1 ) =C ( 1 ) / ( 2*P I ) ELSE C ( I ) =C ( I ) / (2*PI ) +C ( 1-1 ) 

1480 GOSUB 1670 

1490 NEXT I: CL0SE#2: GOTO 1770 

1500 OPEN OUT 1 $ FOR OUTPUT AS 1: OPEN 0UT2$ FOR OUTPUT AS 2 

1510 PR I NT # 1 , "Andrus Stop-Look Detection Algori thm" , DATES: PR INT# 1 ,: PR INT# 1 , 
1520 PR INT# 1 t "NUMBER OF (X,Y) =" ; : PRINTtl , USING"#####" 5 N; 

1530 PRINT#1 , , " I : SEARCH POINT NUMBER" 

1540 PRINTI1 , "TIME LATE =";:PRINT#1, US ING" #####.##"; TL ; 

1550 PR INT#1 , , " X , Y : SEARCH POINT COORDINATES" 

1560 PR I NT# 1 , "SEARCH VELOCITY =";:PRINT#1, USING"#####. ##" 5 SV; 

1570 PRINT#l,, n TIME : TIME SEARCHER AT (X,Y)" 

1580 PR I NT# 1 1 "TARGET VELOCITY =";:PRINT#1, USING "#####. ##" 5 TV; 

1590 PRINT#!,," TLCR : TARGET LOCATION CIRCLE RADIUS" 
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1600 PR INT#1 , "DETECTION RANGE =";:PRINT# 1 , USING"#####. ##” ;R; 

1610 PRINT#1 , , " A1,A2: TARGET/SEARCH CIRCLE INTERSECTION ANGLES" 

1620 PR I N T# 1 , , , " C s CUMULATIVE PROBABILITY OF DETECTION" 

1630 PRINT#1 , : PRINT# 1 ,: PRINT#2, 

1640 PRINT#1 , " I X Y TIME TLCR A1 A2 C" 

1650 PRINT#2 , " B( I )= CUMULATIVE ANGLE DETECTION COVERAGE AT X ( I ) , Y ( I ) . " 
1660 PRINT# 1 , : RETURN 710 
1670 P R I N T# 1 , USING "###";I; 

1680 PR I NT# 1 , USING "######. ##";X(I) ,Y(I) ,TI (I) ,TR(I) ,A1 (I) ,A2 (I) ,C(I) 

1690 IF A 1 ( I ) < >A2 ( I ) AND C ( I ) < >C ( I- 1 ) THEN 1700 ELSE RETURN 

1700 L=1 s PRINT#2, : PRINT#2," B ( " ; I ; " ) = " 5 : IF I<10 THEN PRINT#2, " " ; 

1710 FOR J=1 TO K 

1720 PRINT#2, PRINT#2, USING "#.##"; B1 (J ) 5 : PRINT#2, 

PRINT#2 , USING "#.##"; B2 (J ); : PRINT#2, ") 

1730 IF L=5 AND JOK THEN PRINT#2, 

1740 L=L + 1 : IF L>5 AND JOK THEN PRINT#2, " " 5 : L = 1 

1750 NEXT J 
1760 RETURN 

1770 OPEN 0UT2S FOR INPUT AS 2 
1780 WHILE NOT EOF (2) 

1790 LINE INPUT#2 , A$ 

1800 PRINTftl , A$ 

1810 WEND 

1820 ET=TIMER: IF ET < BT THEN A = 24*60*60-BT+ET ELSE A=ET-BT 
1830 PR I NT# 1 , : PR INT#1 , "Execution Time = PRINT#1, USING "####.##“; A/60; s 
PRINT#1 , " Minutes" 

1840 CL0SE#1: OPEN OUT 1 # FOR INPUT AS 1 
1850 CLS 

1860 FOR 1=1 TO 23 

1870 IF EOF ( 1 ) THEN 1920 

1880 LINE I NPU T # 1 , A$ 

1890 PRINT A* 

1900 NEXT I: LOCATE 25,24: PRINT "PRESS ENTER KEY TO SCROLL SCREEN"; 

1910 A$=INKEY$: IF A$="" THEN 1910 ELSE CLS: GOTO 1860 
1920 CLOSE# 1 : CL0SE#2: KILL 0UT2I: END 
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Appendix B: Computation Procedures: 

The computations used in the program of Appendix A tor the Stop 
and Look algorithm are: 

a. S(i> = SQR(X(i) A 2+Y(i)' N 2) . 

S(i) is the distance trom (0,0) to <X(i),Y(i)). 

b. ST (i ) = SUM C SDR (X(j)— X(j — 1) ) A 2+ ( Y ( j ) — Y ( j — 1 ) ) A 2 3 

tor j =1 to i and ST(1)=0. 

ST(i) is the total distance moved by the searcher trom 
(X < 1 > , Y ( 1 ) ) to (X (i ) , Y (i ) ) . 

c. TI(i) = ST ( i ) /SV+TL. 

T(i) is the time the searcher is at (X(i),Y(i)). SV is 
the search velocity and TL is the time late. 

d. TR (i ) = TI <i>*TV. 

TR(i) is the radius ot the target location circle. TV is 
the target velocity. 

e. r as TI (i >*TV. 

The equation tor the target location circle at time 
TI(i). r is the radial coordinate, 
t. R~2 - r A 2-2*r*S(i >*C0S(A-ATAN(Y(i ) /X (i > ) )+S(i)^2. 

The equation tor the searcher detection circle at time 
TI(i). A is angle coordinate. R is the detection range, 
e. A(i) = (Al(i), A2 ( i ) ) 

=AC0S (R~2— S(i )^2— TR(i ) A 2) / ( -2*S ( i ) *TR ( i ) ) + AT AN ( Y ( i ) /X (i) ) . 
A(i) is the pair ot polar angle coordinates tor the 
intersection ot the target location circle and the 
searcher detection circle at time TI(i). 
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